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A primary goal in continuum strong QCD is to 
develop an intuitive understanding of the spec- 
trum and interactions of hadrons in terms of 
QCD's elementary degrees of freedom. In ad- 
dressing this an efficacious strategy is to employ 
a single framework in calculating a prodigious 
number of observables. With the aim being an 
elucidation of hadronic structure and nonpertur- 
bative aspects of QCD, a focus on the electroweak 
interactions of hadrons is most useful because in 
this case the probe is well understood and the 
features of the hadronic target are unambigu- 
ously under scrutiny. Both confinement and dy- 
namical chiral symmetry breaking (DCSB) play 
key roles in determining these features, and are 
the main themes of this discourse. 

As reviewed in Ref. [1], the last decade has 
seen a modest renaissance in the use of Dyson- 
Schwinger equations (DSEs) [2]: in exploring 
their formal foundation and in their phenomeno- 
logical application. They provide the medium for 
this discussion. 

The DSEs are a nonperturbative means of 
analysing a quantum field theory. Derived from 
a theory's Euclidean space generating functional, 
they are an enumerable infinity of coupled inte- 
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gral equations whose solutions are the n-point 
Schwinger functions [Euchdean Green functions] , 
which are the same matrix elements estimated in 
numerical simulations of lattice-QCD. In theo- 
ries with elementary fermions, the simplest of the 
DSEs is the gap equation, which is basic to study- 
ing dynamical symmetry breaking in systems as 
disparate as ferromagnets, superconductors and 
QCD. The gap equation is a good example be- 
cause it is familiar and has all the properties that 
characterise each DSE: its solution is a 2-point 
function [the fermion propagator] while its kernel 
involves higher n-point functions; e.g., in a gauge 
theory, the kernel is constructed from the gauge- 
boson 2-point function and fermion-gauge-boson 
vertex, a 3-point function; a weak-coupling ex- 
pansion yields all the diagrams of perturbation 
theory; and solved self-consistently, the solution 
of the gap equation exhibits nonperturbative ef- 
fects unobtainable at any finite order in pertur- 
bation theory; e.g, dynamical symmetry break- 
ing. 

The coupling between equations; i.e., the fact 
that the equation for a given m-point function al- 
ways involves at least one n > m-point function, 
necessitates a truncation of the tower of DSEs 
in order to define a tractable problem. One sys- 
tematic and familiar truncation is a weak cou- 
pling expansion to reproduce perturbation the- 
ory. However, that precludes the study of non- 
perturbative phenomena and hence something 
else is needed for the investigation of strongly in- 
teracting systems and bound state phenomena. 

In analysing the ferromagnetic transition, the 
Hartree-Fock approximation yields qualitatively 
reliable information and in QCD its analogue: 
the rainbow truncation, has proven efficacious. 
However, a priori it can be difficult to judge 
whether a given truncation will yield reliable re- 
sults and a systematic improvement is not always 
obvious. It is here that some model-dependence 
enters but that is not new, being typical in the 
study of strongly- interacting few- and many-body 
systems. 

To proceed with the DSEs one just employs 
a truncation and explores its consequences, ap- 
plying it to different systems and constraining 
it, where possible, by comparisons with experi- 
mental data, and also with other theoretical ap- 



proaches on their common domain of application. 
In this way a reliable truncation can be iden- 
tified, and then attention paid to understand- 
ing the keystone of its success and improving its 
foundation. This pragmatic approach has proven 
rewarding in strong QCD, not least because a 
correctly-executed weak coupling expansion is 
guaranteed to match onto perturbation theory 
so that modelling is restricted to the infrared do- 
main. 

2. Quark Propagator 

Obviously it is only possible to study DCSB in 
theories with a well-defined chiral limit. Asymp- 
totically free theories such as QCD are in this 
class. A useful starting point for any discus- 
sion of DCSB is the renormalised quark-DSE [see 
App. ^for the Euclidean metric conventions used 
herein] : 



S{p) ^ = Z2 [i-i ■p + mbaro) 



(2.1) 



+ g^D^,{p-q)^^^S{q)Tl{q,p), 

where D^„{k) is the renormalised dressed-gluon 
propagator, T'^^{q;p) is the renormalised dressed- 
quark-gluon vertex, mbarc is the A-dependent 
current-quark bare mass that appears in the La- 
grangian and := /'^ (i^g/(27r)'* represents 
mnemonically a translationally-invariant regular- 
isation of the integral, with A the regularisation 
mass-scale. The final step in any calculation is 
to remove the regularisation by taking the limit 
A ^ 00. 

Using a translationally invariant regularisa- 
tion makes possible the preservation of Ward- 
Takahashi identities, which is crucial; e.g., in 
studying DCSB [3]. One implementation well- 
suited to a nonperturbative solution of the DSE 
is Pauli-Villars regularisation, in which quarks 
interact with an additional massive gluon-like 
vector boson: mass ~ A, that decouples as A ^ 
00 [4]. An alternative is a numerical implemen- 
tation of dimensional regularisation, which, al- 
though more cumbersome, can provide the nec- 
essary check of scheme-independence [5] . 



In Eq. (|2j), Zi(C^A^) and Z2(C,A^) are 



the quark-gluon- vertex and quark wave function 
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renormalisation constants, which depend on the 
renormahsation point, C,, and the regularisation 
mass-scale, as does the mass renormahsation con- 
stant 

z„,(c^ = Z4(C^A2)/Z2(c^ A2), (2.2) 

with the renormahsed mass given by 

m(C) :=TObare(A)/^m(C',A2). (2.3) 

Although I have suppressed the flavour label, S, 
r° and mbaro depend on it. However, one can al- 
ways use a flavour-independent renormalisation 
scheme, which I assume herein, and hence all 
the renormalisation constants are flavour-inde- 
pendent [4]. 



The solution of Eq. (2.1) has the form 
S{p) ^ 



z{p\e) 

i-fp + M{jP,C) ' 



(2.4) 
(2.5) 



The functions ^(p^, C^), B(p'^, ^^) embody all the 
effects of vector and scalar quark-dressing induced 
by the quark's interaction with its own gluon 
field. The ratio: A/(p^,C^), is the quark mass 
function and a pole mass; i.e., the on-shell mass, 
would be the solution of 



'polo 



M2(p2 



i2„i„C')-0. (2.6) 



A widely posed conjecture is that confinement 
rules out a solution of this equation [6] , and that 
is discussed further in Sec. 4 



Equation (2.1) must be solved subject to a 
renormalisation [boundary] condition, and 
because the theory is asymptotically free it is 
practical and useful to impose the requirement 
that at a large spacelike 



s{py 



17 • p -I- m(C) , 



(2.7) 



where m{C,) is the renormahsed current-quark 
mass at the scale C,. By "large" here I mean 
C'^ S> Aq(^j-, so that in quantitative, model stud- 
ies extensive use can be made of matching with 
the results of perturbation theory. It is the ultra- 
violet stability of QCD; i.e., the fact that pertur- 
bation theory is valid at large spacelikc momenta, 
that makes possible a straightforward definition 



of the chiral limit. It also provides the starkest 
contrast to strong-coupling QED, whose rigorous 
definition remains an instructive challenge [2,7]. 

Multiplicative renormalisability in gauge the- 
ories entails that 



A(p^c^) _ z,(c^A^) 



1 



A(C2,C2) 
(2.8) 

and beginning with Ref. [8] this relation has been 
used efficaciously to build realistic Ansalze for 
the fermion-photon vertex in quenched QED. A 
systematic approach to such nonperturbative im- 
provements is developing [9] and these improve- 
ments continue to provide intuitive guidance in 
QED, where they complement the perturbative 
calculation of the vertex [10]. They are also use- 
ful in exploring model dependence in QCD stud- 
ies [11]. 

At one loop in QCD perturbation theory 



^2(C',A2) 



«(C^) 



-7-f//9i 



(2.9) 



running strong-coupling is 

«(C^) 



\'N t — if, and at this order the 



i/3iln CVA^cD 



(2.10) 



In Landau gauge: ^ = 0, so Z2 = 1 at one loop 
order. This, plus the fact that Landau gauge is a 
fixed point of the renormalisation group 
[Eq. ( |2.22| )], makes it the most useful covariant 
gauge for model studies. It also underlies the 
quantitative accuracy of Landau gauge rainbow 
truncation estimates of the critical coupling in 
strong QED [12]. In a self consistent solution of 
Eq. (2.1), Z2 7^ 1 even in Landau gauge but, at 
large C^, the C-dependence is very weak. How- 
ever, as will become evident, in studies of real- 
istic QCD models this dependence becomes sig- 
nificant for < 1-2 GeV'^, and is driven by the 
same effect that causes DCSB. 

The dressed-quark mass function: M(p^, C^) 
= B{p'^X'^)/Mp'^,0, is independent of the renor- 
malisation point; i.e., with C 7^ C 

M{p^ e) = M{p'X^) := M{p^) , Vp^ (2.11) 

It is a function only of p^/Ag^p, which is an- 
other constraint on models. At one loop order 
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the running [or renormalised] mass 



m(C) = M{e) 



iln 



CVA^ 



QCD 



(2.12) 

7„i — 12/ (33 — 2A^/), where rh is the renormahsa- 
tion point independent current-quark mass, and 



the mass renormahsation constant is, Eq. (2.2) 



a(A2 



(2.13) 



The mass anomalous dimension, 7^, is indepen- 
dent of the gauge parameter to all orders in per- 
turbation theory and for two different quark fla- 
vours the ratio: 



(2.14) 



which is independent of the renormahsation point 
and of the renormahsation scheme. The chiral 
limit is unambiguously defined by 



chiral limit 



0. 



(2.15) 



I reiterate now that a weak coupling expan- 



sion of Eq. (2.1) yields each of the diagrams in 
perturbation theory that contributes to the quark 
self energy. However, every one of those contribu- 
tions to B{p^, C^) is proportional to rh and there- 
fore vanishes in the chiral limit; i.e., B{p^, C^) = 
at every order in perturbation theory. 

One finds, in fact, that in the chiral limit 
there is no scalar mass-like divergence in the cal- 
culation of the self energy. This is manifest in the 



quark DSE, with Eq. (2.1) capable of yielding, in 
addition to the perturbative result: B{p'^,(^) = 
0, a solution M{p^) = B{p^,C^)/A{P^,C^) 7^ 
that is power-law suppressed in the ultraviolet: 
M(jP) ~ guaranteeing convergence of the 

associated integral without subtraction. This is 
dynamical chiral symmetry breaking 

DCSB : AI{p^) ^ when m = . (2.16) 

As we shall see, in QCD this is possible if and 
only if the quark condensate is nonzero: the cri- 
teria are equivalent, and its existence places con- 



straints on the kernel in Eq. (2T), as discussed 
further in Sec. ^. 



2.1 Gluon Propagator 

That kernel is constructed from the dressed-gluon 
propagator and the dressed-quark-gluon 



vertex, and encodes in Eq. (2.1) all effects of the 
quark-quark interaction. In a covariant gauge 
the renormalised dressed-gluon propagator is 

i^M^lKj - I p H i; , 

(2.17) 

where d(fc2,C^) = l/[l+n(fc2, C^)], with n(fc2, C^) 
the renormalised gluon vacuum polarisation for 
which the conventional renormahsation condition 
is 



^(c^c') = o; i.e.,d(c^c') 



(2.18) 



For the dressed-gluon propagator, multiplica- 
tive renormalisability entails 



d{k\e) ^3(c^A2) 



1 



(2.19) 

and at one loop in perturbation theory 



^3(C',A^) = 



C,(A2) 1^^i/fa 



(2.20) 



71 ~ \N} — i(13 — 3^). The gauge parameter 
is also renormahsation point dependent; i.e., the 
renormalised theory has a running gauge param- 
eter. However, because of Becchi-Rouet-Stora 
[BRST or gauge] invariance, there is no new dy- 
namical information in that: its evolution is com- 
pletely determined by the gluon wave function 
renormahsation constant 

C(C')=^3"'(C',A')ebare(A). (2.21) 

One can express ^(C^) in terms of a renormahsa- 
tion point invariant gauge parameter: ^, which 
is an overall multiplicative factor in the formula 
and hence 

Landau Gauge : | = ^ ^(C^) = (2.22) 

at all orders in perturbation theory; i.e.. Lan- 
dau gauge is a fixed point of the renormahsation 
group. 
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2.2 Quark-gluon Vertex 

The other element of the kernel, the renormalised 
dressed-quark-gluon vertex, has the form 

Tl{k,p) = ^TAk,p); (2.23) 

i.e., the colour matrix structure factorises. It is a 
fully amputated vertex, which means all the an- 
alytic structure associated with elementary ex- 
citations has been eliminated. To discuss this 
further I introduce the notion of a particle-like 
singularity. 

A Particle-like Singularity is one of the form: 



1 



(p2 + ^2)c 



a e (0,1] 



(2.24) 



If the vertex possesses such a singularity then 
its P-dependence can be expressed via a non- 
negative spectral density, which is impossible if 
a > 1. a = 1 is the ideal case of an isolated 5- 
function distribution in the spectral densities and 
hence an isolated free-particle pole, a £ (0, 1) 
corresponds to an accumulation at the particle 
pole of branch points associated with multiparti- 
cle production, as occurs with the electron prop- 
agator in QED because of photon dressing. 

The dressed-quark-gluon vertex is a fully am- 
putated 3-point function. Therefore in this case 
the presence of such a singularity would entail the 
existence of a flavour singlet composite [quark- 
antiquark bound state] with colour octet quan- 
tum numbers and mass ni — b. [The bound 
state amplitude follows immediately from the as- 
sociated homogeneous Bethe-Salpeter equation 
(BSE), which the singularity allows one to de- 
rive.] However, an excitation like this must not 
exist as an asymptotic state: that would vio- 
late the observational evidence of confinement. 
Hence, as discussed further on page I con- 
clude that the vertex should not exhibit a particle- 
like singularity, and any modelling of r^(A:,p) 
ought to be consistent with this constraint. [NB. 
a > 1 yields an admissible non-particle-like sin- 
gularity.] 

Expressing the Dirac structure ofT^{k,p) re- 
quires twelve independent scalar functions: 



which it is not necessary to articulate fully herein. 
A pedagogical discussion of the perturbative cal- 
culation of r,y{k,p) can be found in Ref. [13] 
while Refs. [14,15] explore its nonperturbative 
structure and properties. I only make Fi{k,p,Q 
explicit because the renormalisability of QCD en- 
tails that it alone is ultraviolet divergent. Defin- 
ing 

/i(fc2,C') :=Fi(A:,-A:,C), (2.26) 

the conventional rcnormalisation boundary con- 
dition is 

/i(C',C') = l, (2.27) 

which is practical because QCD is asymptotically 
free. Multiplicative renormalisability entails 



/l(fc^c^) 
Zl(c^A2) 



(2.28) 



Zi(C2,A2) 
and at one loop order 



/l(c^c^) 



-1 -7r//3i 



(2.29) 



Zi(/x2,A2) 

7^ = i[!(3-l-0^ 
3. DCSB 

At this point each element in the quark DSE, 



r,(fc,p) = 7,Fi(fc,p,c) + 



(2.25) 



Eq. (2.1), is defined, with some of their pertur- 
bative properties elucidated, and the question 
is: "How does that provide an understanding of 
DCSB?" It is best answered using an example, in 
which the model-independent aspects are made 
clear. 

The quark DSE is an integral equation and 
hence its elements must be known at all val- 
ues of their momentum arguments, not just in 
the perturbative domain but also in the infrared. 
While the gluon propagator and quark-gluon ver- 
tex each satisfy their own DSE, that couples the 
quark DSE to other members of the tower of 
equations and hinders rather than helps in solv- 
ing the gap equation. Therefore, as with all ap- 
plications of the gap equation, one employs An- 
satze for the interaction elements [D^v{k) and 
r,y(fc,p)], constrained as much and on as large 
a domain as possible. This approach has a long 
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history in exploring QCD [2] and I illustrate it 
using the model of Ref . [4] . 

The renormalised dressed-ladder truncation 
of the quark-antiquark scattering kernel [4-point 
function] is 



Tl{p+,q+)S{q+) S{q^)Tl{q^,p^) 



(3.1) 



where p± = p ± P/2, q± ^ q ± P/2, with P 
the total momentum of the quark-antiquark pair, 
and although I use it now I have suppressed the 
^-dependence of the Schwingcr functions. From 
Eqs. (|2.18H2.20| ) it follows that for [p-qf 



large and spacelike 



^3(c^A^) 

^3(Q2,A2) 



71 //3i 



(3.2) 
(3.3) 



D';::%P~q). (3.4) 



Using this and analogous results for the other 
Schwinger functions then on the kinematic do- 
main for which ^ p^ ^ q^ is large and space- 
like[g2(C2) :=4™(C2)] 



if(p,g;P)-«4^a(Q2)i^ftec(p_^) 



(3.5) 



because Eqs. (|2j), ( p^ and ( ^^ yield 



/3i Pi 



f3i 



= 1 



(3.6) 



This is one way of understanding the origin 
of an often used Ansatz in studies of the gap 
equation; i.e., making the replacement 



5'^^.(fc)->47ra(fc2)i5fr-(A:) 



(3.7) 



in Eq. {2A), and using the "rainbow truncation:'' 



(3.8) 



Equation (3.7) is often described as the "Abe- 
lian approximation" because the left- and right- 
hand-sides [r.h.s.] are equal in QED. In QCD, 
equality between the two sides cannot be ob- 
tained easily by a selective resummation of di- 
agrams. As reviewed in Ref. [2], Eqs. (5.1-5.8), 



it can only be achieved by enforcing equality be- 
tween the renormalisation constants for the 
ghost-gluon vertex and ghost wave function: Zi 
= Z3. A mutually consistent constraint, which 
follows formally from Zi — Z3, is to enforce 
the Abelian Ward identity: Zi — Z2. At one- 
loop this corresponds to neglecting the contribu- 
tion of the 3-gluon vertex to T^, in which case 
7r ^ §5 = jF- This additional constraint pro- 
vides the basis for extensions of Eq. (3.8); i.e.. 



using Ansdtze for that are consistent with the 
QED vector Ward-Takahashi identity; e.g., Refs. 
[11]. 

Arguments such as these inspire the follow- 
ing Ansatz for the kernel in Eq. ( |2.l|) [4] : 

Zi g'D^,ip-q)^j,,S{q)rtiq,p) (3.9) 

- J^^ Qiip - qf) D'^Tip - 9) Y7M5(g) , 

with the ultraviolet behaviour of the so-called 
"effective coupling:" ^?(fc^), fixed by that of the 
running strong-coupling. Since it is not possible 
to calculate Zi nonperturbatively without analy- 
sing the DSE for the dressed-quark-gluon vertex, 
this Ansatz absorbs it in the model effective cou- 
pling. 

Equation ( |3.9| ) is a model for the product of 
the dressed-propagator and dressed-vertex and 
its definition is complete once the behaviour of 
G{k^) in the infrared is specified; i.e., for fc^ < 1- 
2GeV^. Reference [4] used 



In 



with J^(fc2) 

„2 



(1 + F/A 



2 

QCD 



■T{k^) ,{'S.10) 



[1 - exp(-fcV[4TO2])]/fc2 ^ 



1. For Nf = 4, Aq^d* = 0.234 GeV. 



The qualitative features of Eq. ( 3.1C| ) are 
plain. The first term is an integrable infrared 
singularity [16] and the second is a finite-width 
approximation to S'^{k), normalised such that it 
has the same / d'^k as the first term. In this 
way the infrared strength is split into the sum of 
a zero-width and a finite-width piece. The last 



term in Eq. ( 3.10 ) is proportional to a{k )/k at 
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large spacelike-A:^ and has no singularity on the 
real-A:^ axis. 

There are ostensibly three parameters in 



When m 7^ the renormalisation condition, 



Eq. (3.10): D, u) and mt- However, in Ref. [4] 
uj = 0.3GeV(= l/[.66fm]) and m* = 0.5GeV(= 
l/[.39 Im]) were fixed, and only D and the renor- 
malised u = d- and s-current-quark masses were 
varied in an attempt to obtain a good description 
of low-energy tt- and X-meson properties, using 
a renormalisation point C = 19 GeV that is large 
enough to be in the perturbative domain. [The 
numerical values of lo and mt are chosen so as to 
ensure that ^(fc^) sa 47ra(fc2) for fc^ > 2Ge\f'^. 
Minor variations in uj and mt can be compen- 
sated by small changes in D.] Such a proce- 
dure could self-consistently yield Z) = 0, which 
would indicate that agreement with observable 
phenomena precludes an infrared enhancement 
in the effective interaction. However, that was 
not the case and a good fit required 



D = (0.884 GeV)2, 
with renormalised current-quark masses 



(3.11) 



m„,d(C) = 3.74 MeV, m,(C) = 82.5 MeV, 

(3.12) 

which are in the ratio 1 : 22, and yielded, in MeV, 







mx 




Ik 


Calc. 


[4] 


139 


497 


131 


154 


Expt. 


[17] 


139 


496 


131 


160 



(3.13) 



and other quantities to be described below. An 
explanation of how this fit was accomplished re- 
quires a discussion of the homogeneous Bethe- 
Salpeter equation, which I postpone until Sec. 0. 
It is described in detail in Refs. [4,18]. Here I 
focus instead on describing the properties of the 
DSE solution obtained with these parameter val- 
ues. 

Using Eqs. (2.1-2!3|) and ( |3.9| ) the gap equa- 
tion can be written 

S(pXr^ ^ Z2ii-p+Zim{C) + Y.'{p,K) , (3.14) 

with the regularised quark self energy 

I]'(p,A):= (3.15) 



Eq. ( |2.7|) , is straightforward to implement. Writ- 
ing 

S'(p, A) *7 .p {A'{p\K^) - l)+B'{p\K^) , 

(3.16) 

which emphasises that these functions depend on 
the regularisation mass-scale. A, Eq. ( ^.7| ) entails 



^2(C',A2) 



A'(C^A2), 



(3.17) 



m(C) = Z2(C^ A^) mb„.(A2) -f B'(C^ A^) (3.18) 
so that 

A{p\e) = l + A'{p\K^) - A'(C^ A^) , (3.19) 

B{p\e) = m(C) + B'{p\K^) ~ i?'(C^ A2) . 

(3.20) 

Multiplicative renormalisability requires that 
having fixed the solutions at a single renormali- 
sation point, their form at another point, C, is 
given by 

S-\pX) = (3.21) 

This feature is evident in the solutions obtained 
in Ref. [4]. It means that, in evolving the renor- 



malisation point to C, the "1" in Eqs. ( 3.19 ) is re- 
placed by Z2(C^ A2)/Z2(C^ A2), and the "m(C)" 
by m((); i.e., the "seeds" in the integral equa- 
tion evolve according to the QCD renormalisa- 



tion group. This is why Eq. ( p.9| ) is called a 
"renormalisation-group-improved rainbow trun- 
cation." 

Returning to the chiral limit, it follows from 
Eqs. (Ill), (|3|), ( pl| ) and ( |2l5| ) that for m = 



^2(C',A^)mbarc(A2)=0, VA. 



(3.22) 



Hence, as remarked on page there is no sub- 
traction in the equation for B{p^,(^): 



i.e., Eq. (|3.19| ) becomes 



B{p',C ) = B'{p^, A") , lim B'{p\K') < oo , 

(3.23) 

which is only possible if the mass function is at 
least 1/p^-suppressed. This is not the case in 
quenched strong-coupling QED, where the mass 
function behaves as [19] 



cx cos(const. ln[p''/C ])/(p /C 



a2/^2Nl/2 



(3.24) 
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Figure 1: Quark mass function obtained as a solu- 



tion ofEq. (2.1) using the model ofEqs. (3.9), (3. IC) 



and current-quark masses, fixed at a renormalisation 

= 82 MeV, 
The indi- 



point C = 19GeV: m^^ = 3.7 MeV, 
mi = 0.58 GeV and = 3.8 GeV 
cated solutions of M^{p^) = define the Euclidean 
constituent-quark mass, Mf , which takes the values: 
= 0.56 GeV, Mf = 0.70 GeV, Mf = f.SGeV, 
Mjf — 4.6 GeV. These values and their ratios are 
consistent with contemporary phenomenology; e.g., 
Refs.[f7,20] (Figure adapted from Refs. [2f ,22].) 



and that is the origin of the complications dis- 
cussed in Refs. [2,5,7]. 

In Fig. |l| I present the renormalised dressed- 
quark mass function, M(p^), obtained by solving 
Eq. ( |3.14[ ) using the model and parameter values 
of Ref. [4], Eqs. (|j[f]T|), and also in the chiral 
limit and with typical heavy-quark current-mass 
values. 

In the presence of explicit chiral symmetry 
breaking Eq. ( 2.12 ) describes the form of M(p^) 



for > 0(1 GeV^). In the chiral limit, however, 
the ultraviolet behaviour is given by 



2^2^„ 



(3.25) 



[\ In 



QCD 



1-7^ 



where {qq)'^ is the renormalisation-point-indepen- 
dent vacuum quark condensate. This behaviour 
too is characteristic of the QCD renormalisation 
group [23] and exhibits the power-law suppres- 
sion anticipated on page ^. These results for 
the large-p2 behaviour of the mass function are 



model independent; i.e., they arise only because 
the DSE truncation is consistent with the QCD 
renormalisation group at one loop. (It has long 
been known that the truncation defined by 
Eq. ( |3.9| ) yields results in agreement with the 
QCD renormahsation group at one loop; e.g., 
Refs. [24].) 

The gauge invariant expression for the renor- 
malisation-point-dependent vacuum quark con- 
densate was derived in Ref. [25]: 

r-A 



Z4(C^A2)7VetrD / 5°(g,C), 

J a 



(3.26) 

where tr^i identifies a trace over Dirac indices 
only and the superscript "0" indicates the quan- 
tity was calculated in the chiral limit. Substi- 



tuting Eq. ( |3.25D into Eq. ( |3.26| ), recaUing that 
Zi = Zra in Landau gauge and using Eq. ( 2.13| ) 
leads to the one-loop expression 



(3.27) 



Employing Eq. ( |2.12|) , this exemplifies the gen- 
eral result that 



(3.28) 



i.e., that this product is renormalisation point 
invariant and, importantly, it shows that the be- 



haviour expressed in Eq. ( p. 251) is exactly that 
required for consistency with the gauge invariant 
expression for the quark condensate. A model, 
such as Ref. [26], in which the scalar projection 
of the chiral limit dressed-quark propagator falls 
faster than l/p"*, up to In-corrections, is only con- 
sistent with this quark condensate vanishing, and 
it is this condensate that appears in the current 
algebra expression for the pion mass [25], as dis- 
cussed in connection with Eq. (7 



Equation (3.25) provides a reliable means of 
calculating the quark condensate because cor- 
rections are suppressed by powers of Aqq^/^^. 
Analysing the asymptotic form of the numerical 
solution one finds 



-(w)" = (0.227 GeV)' 



(3.29) 



Using Eq. ( |3.27| ) one can define a one- loop evolved 
condensate 

-(9'?>clc=lGcV ■= (1^^ [1/Aqcd])^" {qq? , 

(3.30) 
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which takes the value 



GcV 



(0.241 GeV)= 



(3.31) 



This can be directly compared with the value of 
the quark condensate employed in contemporary 
phenomenological studies [27]: 



(0.236 ± 0.008 GeV)^ 



(3.32) 



It was noted in Ref. [4] that increasing uj 
1.5 w in G{k'^) increases the calculated value in 
Eq. ( 3.30 ) by ^ 10%; i.e., the magnitude of the 
condensate is correlated with the degree of in- 
frared enhancement/strength in the effective in- 
teraction. That is unsurprising because it has 
long been known that there is a critical cou- 
phng for DCSB; i.e., the kernel in the gap equa- 
tion must have an integrated strength that ex- 
ceeds some critical value [24]. This is true in all 
fermion-based studies of DCSB, a point discussed 
in more detail on page ^ 

The renormalisation-point-invariant current- 
quark masses corresponding to the mf{Q in 
Fig. ^ ar e obtained in the following way: using 
Eq. (3.2(;), direct calculation from the chiral limit 
numerical solution gives 



C=19 GoV 



(0.275 GeV)^ 



(3.33) 



and hence from the values of = mf{C,) listed 



in Fig. g and Eqs. ( |3.28D , ( p^ ), in MeV 



rhu.d = 6.60 , rhs = 147 , 
rhc — 1 030 , "nib — 6 760 , 



(3.34) 



from which also follow one-loop evolved values in 
analogy with Eq. ( |3.30|) : 



„lGcV _ 
,lGcV 



5.5, ml^''^ = 130, 
860, ml 



IGcV 



5 700. 



(3.35) 



Smaller values of the light-current-quark masses 
would require a larger value of the vacuum quark 
condensate in order to be consistent with light- 
meson properties. 

Figure ^ highlights a number of qualitative 
aspects of the quark mass function. One is the 
difference in the ultraviolet between the behavi- 
our of M(p^) in the chiral limit and in the pres- 
ence of explicit chiral symmetry breaking. In the 
infrared, however, the m, d-quark mass function 



and the chiral limit solution are almost indistin- 
guishable. The chiral limit solution is nonzero 
onlyhecsLUse of the nonperturbative DCSB mech- 
anism whereas the u, d-quark mass function is 
purely perturbative at > 20 GeV^. Hence the 
evolution to coincidence of the chiral-limit and 
u, d-quark mass functions makes clear the tran- 
sition from the perturbative to the nonperturba- 
tive domain. 

It is on this nonperturbative domain that 
A(p^,C^) differs significantly from one. A con- 
comitant observation is that the DCSB mecha- 
nism has a significant effect on the propagation 
characteristics of u, d, s-quarks. These aspects 
of the momentum-dependence of the scalar func- 
tions in the dressed-light-quark propagators have 
recently been confirmed in numerical simulations 
of lattice-QCD [28], as illustrated in the compar- 
isons of Fig. H 

Returning to Fig. |^, it is evident that the 
large current-quark mass of the 6-quark almost 
entirely suppresses momentum-dependent dress- 
ing, so that A/f,(p^) is nearly constant on a sub- 
stantial domain. This is true to a lesser extent 
for the c-quark but in both cases one observes 
that DCSB is not the dominant mass-generating 
mechanism. 

A means of quantifying the effect of the 
DCSB mechanism on dressed-quark propagation 
characteristics is now obvious. One can introduce 
[21,22]: Cf := Mf/m'), where Mf is the Eu- 
clidean constituent-quark mass, defined in Fig. |l| 
as the solution of 

{Mff - IvP{p^ = (Aff )^ C^) = , (3.36) 

and this ratio takes the values 



/ 



-/ 



150 10 2.2 1.2 



(3.37) 



The values are representative and definitive: 



for light-quarks Lq=u,d,. 



10-100, while for 



heavy-quarks >Cq=c,6 ''^ 1; and highlight the exis- 
tence of a mass-scale characteristic of DCSB: M^. 
The propagation characteristics of a flavour with 
rrij < are significantly altered by the DCSB 
mechanism, while for fiavours with 3> 
momentum-dependent dressing is almost irrele- 
vant. It is apparent and unsurprising that ~ 
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ma=0.0603 




0.0 



0.5 



1.0 
pa 
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Figure 2: Upper Panel: Lattice result for Z{pa), 
a ~ 2.0 GeV~^ is tlie lattice spacing, calculated with 
ma = 0.0603 [28], compared with the DSE analogue: 
Z{x), j: is a dimensionless momentum variable, from 
the phenomenological algebraic model of Ref. [29]. 
The variation between the lattice results is an indica- 
tion of extant numerical artefacts. The "dip" in Z{x) 
for a; ~ is a long-standing prediction of phenomeno- 
logically efficacious DSE studies. Lower Panel: Anal- 
ogous plot for the mass function. Numerical arte- 
facts are more significant in this case. Nevertheless 
the enhancement at small-pa, indicative of DCSB, is 
clearly evident and in semi-quantitative agreement 
with that predicted by DSE studies. The dashed-line 
assists with estimating the Euclidean constituent- 
quark mass; ~ 200 MeV in the lattice simulation and 
~ 350 in the phenomenological DSE model [29]. 



0.2 GeV^ Aqcd- This feature of the dressed- 
quark mass function provides the foundation for 
a constituent-quark-hke approximation in the 



treatment of heavy-meson decays and transition 
form factors [22,30]. 

4. Confinement 

To proceed it is necessary to discuss confinement 
and I begin with a definition. Confinement is 
the failure to directly observe coloured excita- 
tions in a detector: neither quarks nor gluons nor 
coloured composites. The contemporary hypoth- 
esis is stronger; i.e., coloured excitations cannot 
propagate to a detector. To ensure this it is suffi- 
cient that coloured n-point functions violate the 
axiom of reflection positivity [31], which is guar- 
anteed if the Fourier transform of the momentum- 
space n-point Schwinger function is not a positive- 
definite function of its arguments. 

Reflection positivity is one of a set of five ax- 
ioms that must be satisfied if the given n-point 
function is to have a continuation to Minkowski 
space and hence an association with a physical, 
observable state. If an Hamiltonian exists for 
the theory but a given n-point function violates 
reflection positivity then the space of observable 
states, which is spanned by the eigenstates of the 
Hamiltonian, does not contain anything corre- 
sponding to the excitation(s) described by that 
Schwinger function. 

The free boson propagator does not violate 
reflection positivity: 



A(x) := 



(2^ 



1 



47r2; 



■ Ki{mx) . 



(4.1) 
(4.2) 



Here x \= {x ■ xY^"^ > and Ki is the monotoni- 
cally decreasing, strictly convex-up, non-negative 
modified Bessel function of the second kind. The 
same is true of the free fermion propagator: 



S{x) 



m 



(2^^ 
Ki{mx) 



Ak-x 



Z7 • k 



7 • X 



(4.3) 



K2{mx) , (4.4) 



which is also positive definite. 

The spatially averaged Schwinger function is 
a particularly insightful tool [32,33]. Consider 
the fermion and let T — x^ represent Euclidean 
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"time," then 



di 



cos{£T) 



-mT 



(4.5) 



Hence the free fermion's mass can easily be ob- 
tained from the large T behaviour of the spatial 
average: 



mT = — lim ln(Ts(r) 



(4.6) 



[The boson analogy is obvious.] This is just the 
approach used to determine bound state masses 
in simulations of lattice-QCD. 

For contrast, consider the dressed-gluon 2- 
point function obtained with [34] 



d(fc2) 



in Eq. ( pTzD 

D{x) : 



Ak-x 



(27r)4 " k^+-f 
dz 



(4.7) 

(4.8) 
(4.9) 



z—'yx 



where ker(z) is the oscillatory Thomson function. 
D{x) is not positive definite and hence a dressed- 
gluon 2-point function that vanishes at A:^ = 
violates the axiom of reflection positivity and is 
therefore not observable; i.e., the excitation it 
describes is confined. 

At asymptotically large Euclidean distances 



cosi^ix + f ) + sin(-i57a; + f ) 



(4.10) 



Comparing this with Eq. (4.1) one identifies a 



mass as the coefficient in the exponential: 

mD=l/V2- (4.11) 

[NB. At large x, Ki{x) cx exp{—x)/^x.] By an 
obvious analogy, the coefficient in the oscillatory 
term is the lifetime [34]: 



l/uiD- 



(4.12) 



Both the mass and lifetime are tied to the dy- 
namically generated mass-scale 7, which, using 

z ,1 ,1 



^2 _|_ ^4 '-^ z + 47^ ^ Z — i"f^ 



(4.13) 



is just the displacement of the complex conjugate 
poles from the leal-k^ axis. 

It is a general result that the Fourier trans- 
form of a real function with complex conjugate 
poles is not positive definite. Hence the existence 
of such poles in a n-point Schwinger function is a 
sufficient condition for the violation of reflection 
positivity and thus for confinement. 

The spatially averaged Schwinger function is 
also useful here. 



D{T) := J d^xD{x,T) 



= ^e-727T^os(^7T+f), (4.14) 



and, generalising Eq. (4.6), one can define a T- 
dependent mass: 



m 



i{T)T := -\nD{T) 
= ln(27) + -^7r-ln 



cos(^7r+f 



(4.15) 



It exhibits periodic singularities whose frequency 
is proportional to the dynamical mass-scale that 
is responsible for the violation of refiection posi- 
tivity. If a dressed-fermion 2-point function has 
complex conjugate poles it too will be charac- 
terised by a T-dependent mass that exhibits such 
behaviour. 

4.1 Three-dimensional QED 

This reflection positivity criterion has been em- 
ployed to very good effect in three dimensional 
QED [35]. First, some background. QED3 is 
confining in the quenched truncation [36]. That 
is evident in the classical potential 

d^k 



Vir) :- 



dx:i 



(27r)a 



^ik-x+iksxs 



= ^ln(eM, r^=xl+4, (4-16) 

which describes the interaction between two static 
sources. [NB. has the dimensions of mass in 
QED3.] It is a logarithmically growing poten- 
tial, showing that the energy required to sepa- 
rate two charges is infinite. Furthermore, V{r) 
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is just a one-dimensional average of the spatial 
gauge-boson 2-point Schwinger function and it 
is not positive definite, which indicates that the 
photon is also confined. 

However, if now the photon vacuum polari- 
sation tensor is evaluated at order using Nf 
massless fermions then, with the notation of 



Eq. ( 2.17 ), the photon propagator is characterised 
by [37] 



d(fc2) 



1 



from n(fc^ 



a 
k 



(4.17) 



fc2 + ak 

a = iV/e^/S, and one finds [38] 

Vir)^^[Uoi&r)~Noi&r)] , (4.18) 

where Ho(a;) is a Struve function and No{x) a 
Neumann function, both of which are related to 
Bessel functions. In this case V{r) is positive 
definite, with the limiting cases 



V{r) - ln(<5r) , V{r) 



TT^^ (4-19) 



and confinement is lost in QED3. That is easy 
to understand: pairs of massless fermions cost 
no energy to produce and can propagate to in- 
finity so they are very effective at screening the 
interaction. 

With d(fc2) = 1/[1 + n(fc2)] and sensible, 
physical constraints on the form of n(fc^), such 
as boundedness and vanishing in the ultraviolet, 
one can show that [39] 



1 



2.i + n(0)^"(^'^) + ^°"^*- + "(^)' 

(4.20) 

where h{r) falls-off at least as quickly as l/r. 
Hence, the existence of a confining potential in 
QED3 just depends on the value of the vacuum 
polarisation at the origin. In the quenched trun- 
cation, H(0) — and the theory is logarithmi- 
cally confining. With massless fermions, 
1/[1 -f- H(0)] — and confinement is absent. Fi- 
nally, when the vacuum polarisation is evaluated 
from a loop of massive fermions, whether that 
mass is obtained dynamically via the gap equa- 
tion or simply introduced as an external parame- 
ter, one obtains H(0) < 00 and hence a confining 
theory. 







quenched QED3 o 




massless fermions in vac. pel. + 




fermions with fixed mass in vac.pol. x 




dynamical fermions in vac.pol. □ 













+ 




□ 



Figure 3: A(r) — -crs(T) from Eq. for 
QED3 with 2 flavours of fermion. "Diamond:" con- 
fining, quenched theory; "plus," massless fermions 
used to evaluate the photon vacuum polarisation ten- 
sor; "cross," as before but with fixed-mass fermions; 
"box," same again but fermions with a momentum- 
dependent mass function. (Adapted from Ref. [35].) 



In Ref. [35] the QED3 gap equation is solved 
for all four cases and the fermion propagator anal- 
ysed. The results are summarised by Fig. 0. 



In the quenched theory, Eq. ( 4.16 ), the dressed- 
fermion 2-point function exhibits exactly those 
periodic singularities that, via Eq. dil^ ), are m- 
dicative of complex conjugate poles. Hence this 
feature of the 2-point function, tied to the vi- 
olation of reflection positivity, is a clear signal 
of confinement in the theory. That is empha- 
sised further by a comparison with the theory 
that is unquenched via massless fermions in the 
vacuum polarisation, Eq. (4.17). As I have de- 



scribed, that theory is not confining and in this 
case (Js{T) has the noninteracting, unconfincd 
free particle form in Eq. (4.5). The difference 
could not be more stark. The remaining two 
cases exhibit the periodic singularities that sig- 
nal confinement, just as they should based on 



Eq. ( |4.2q ). 

At this point I note that any concern that 
the presence of complex conjugate singularities 
in coloured n-point functions leads to a violation 
of causality is misguided. Microscopic causality 
only constrains the commutativity of operators, 
and products thereof, that represent elements in 
the space of observable particle states; i.e., the 
space spanned by eigenstates of the Hamiltonian. 
Since Schwinger functions that violate reflection 
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positivity do not have a continuation into that 
space there can be no question of violating causal- 
ity. It is only required that 5-matrix elements 
that describe colour-singlet to colour-singlet tran- 
sitions should satisfy the axioms, including reflec- 
tion positivity. 

The violation of reflection positivity by 
coloured n-point functions is a sufficient condi- 
tion for confinement. However, it is not neces- 
sary, as the example of planar, two-dimensional 
QCD shows [40]. There the fermion two-point 
function exhibits particle-like singularities but 
the colour singlet meson bound state amplitudes, 
obtained from a Bethe-Salpeter equation, van- 
ish at momenta coincident with the constituent- 
fermion mass shell. This excludes the pinch sin- 
gularities that would otherwise lead to bound 
state break-up and liberation of the constituents. 
It is a realisation of confinement via a failure of 
the cluster decomposition property [CDP] [31,41]. 

The CDP is a requirement that the differ- 
ence between the vacuum expectation value of a 
product of fields and all products of vacuum ex- 
pectation values of subsets of these fields must 
vanish faster than any power. [This is modified 
slightly in theories, like QED, with a massless, 
asymptotic state: the photon in this case.] It 
can be understood as a statement about charge 
screening and its failure means that, irrespective 
of the separation between sources, the interaction 
between them is never negligible. That is an ap- 
pealing, intuitive representation of confinement. 
Failure of the CDP is an implicit basis for con- 
finement in the bulk of QCD potential models; 
e.g., Refs. [42]. 

5. Gap Equation's Kernel 

Strong interaction phenomena are characterised 
by DCSB and colour confinement. At low en- 
ergy, DCSB is the more important; for exam- 
ple, in its absence the tt- and p-mesons would be 
nearly degenerate and at the simplest observa- 
tional level that would lead to a markedly differ- 
ent line of nuclear stability. These phenomena 
can be related to the infrared behaviour of ele- 
mentary Schwinger functions in QCD and in this 
section I elucidate some constraints they place 
on this behaviour. 



In Eq. ( |3.9| ) I described an Ansatz for the 
kernel in the quark DSE and used it to elabo- 
rate on the phenomenon of DCSB, arguing that 
a good description of light-meson observables re- 
quires the kernel to exhibit a significant enhance- 
ment in the infrared, Eq. ( ^.11 ) . An obvious ques- 
tion is: "How far-reaching is this result?" 



In general, as is clear from Eq. (3^), the ker- 
nel is a product of two terms: the dressed-gluon 
propagator and the dressed-quark-gluon vertex. 
For > 1-2 GeV^ a perturbative analysis for 



this product is reliable and Eq. (3.9) becomes 
an identity with ^(fc^) — > 47ra(A:^). This means 
that any model-dependence in the Ansatz is con- 
strained to the infrared domain: fc^ < 1-2 GeV^. 

In contemporary DSE applications to QCD 
it is common to build Ansdtze for the higher 
n-point dressed-Schwinger functions and employ 
them in developing intuition about the simpler 
functions. Section |3| provides an illustration. In 
pursuing this certain constraints must be obeyed. 

Of particular interest here is the dressed- 
quark-gluon vertex, which is a fully-amputated 
3-point Schwinger function. It satisfies an inte- 
gral equation that takes the form of an inhomo- 
geneous Bethe-Salpeter equation. The kernel in- 
volves K: the fully-amputated, quark-antiquark 
scattering kernel, which by definition does not 
contain quark-antiquark to single gauge-boson 
annihilation diagrams, such as would describe 
the leptonic decay of the pion, nor diagrams that 
become disconnected by cutting one quark and 
one antiquark line. It also involves the scat- 
tering kernels for: q-q to 2-gluon, K^^; q-q to 
ghost-antighost, K^'^^'^; and q-q to 3-gluon, K^^ , 
and also by definition none of these can contain 
single-gluon intermediate states. Hence, just as 
in the chiral limit a massless pole in the pseu- 
dovector vertex signals the presence of a mass- 
less [pion] bound state, a massless, particle-like 



singularity [see Eq. ( |2.24| )] in T^{q,p) signals the 
presence of a colour-octet bound state in the scat- 
tering matrices: M := K/[l - {SS)K]; Apa 
X29/[i - {DD)K'^3-^; etc. As no such coloured 
bound states has been observed, one must be 
sceptical of calculations or Ansdtze for any of 
the Schwinger functions that entail a particle- 
like singularity in this vertex. [NB. It is inter- 
nally inconsistent to interpret as confined a gluon 



13 



Confinement @ ESI May-July 2000 



CD. Roberts 



whose 2-poiiit function violates reflection posi- 
tivity whilst simultaneously asserting that a par- 
ticle-like singularity in a coloured irreducible 3- 
point function does not describe an asymptotic 
state.] 

The same objection applies to particle-like 
singularities in the fully-amputated, dressed-3- 
gluon vertex, and all like coloured n-point func- 
tions. This anticipates the result of an estimate 
[43] of the 3-gluon vertex via a numerical simu- 
lation of lattice-QCD, which shows no evidence 
for a singularity of any kind. 

Rejecting particle-like singularities in 
r^(g,p), the possibility of an enhancement in the 
kernel of the gap equation can be discussed solely 
in terms of the behaviour of the dressed-gluon 
propagator, which in Landau gauge can be writ- 



ten [cf. Eq. (2.17)] 



D,Ak) = (<5,.-^)A(A:2), (5.1) 



(5.2) 



The question I posed at the beginning of this 
section can now be rephrased as: "Do observable 
strong interaction phenomena necessarily require 

Vik^) > 1 for 1 < fcVAqcD < 10 ?" (5.3) 

[^QCD^ = 234 MeV.] NB. Herein I do not canvass 
the possibility that an irreducible vertex has a 
non-particle-like singularity; i.e., a singularity of 
the form (fc^)~", a > 1. While that evades the 
constraint I have elucidated, there is no indica- 
tion of such behaviour in any study to date. 

I do not have an answer to the question in 
Eq. ( |5.3| ) but the alternatives can be explored. 
The antithesis is the extreme possibility that 



p(fc2=o) = 0, V{k^) <lVfc2 



(5.4) 



which was canvassed in Refs. [34]. ["Extreme" 
because it corresponds to a screening of the fer- 
mion-fermion interaction, as familiar in an elec- 
trodynamical plasma, rather than the antiscreen- 
ing often discussed in zero-temperature chromo- 
dynamics.] 

References [34] proposed solving the DSEs 
via rational polynomial Ansdtze for the one-par- 
ticle irreducible components of the Schwinger 



-1 -I ^ ^ % 

^OT6»OT(p = ''mmms^ - 1/2 Tsl^sMifrrr^ - 1/2 ■-f(r6> Crf«^ 



■ 1/6 '-e5ip-me6^me(\sir- - 1/2 '-trrs^sj-^ "Cx^s^ 



Figure 4: From top to bottom, depictions of the 
DSEs for the gluon (spring), ghost (dashed-line) and 
quark (solid-line) 2-point functions. Following con- 
vention, a filled circle denotes a fully dressed prop- 
agator and an open circle, a one-particle irreducible 
vertex; e.g., the open circle in the first line represents 
the dressed-three-gluon vertex. The figure illustrates 
the interrelation between elements in the tower of 
DSEs: the gluon propagator appears in the DSE for 
the quark and ghost propagator; the ghost and quark 
propagator in the DSE for the gluon, etc. (Adapted 
from Ref. [44].) 



functions. This method attempts to preserve as- 
pects of the organising principle of perturbation 
theory in truncating the DSEs. In connection 
with the gluon DSE, depicted in Fig. ^, these 
studies employed a truncation that, for simplic- 
ity: retains only the first, third and sixth dia- 
grams on the r.h.s. of the first equation in Fig. ^; 
neglects the last [fermion] equation; and employs 
the leading order solution of the ghost equation, 
which has the appearance of the massless free 
propagator: ~ 1/fc^. In this approach the 
Ansdtze for the 3-gluon and quark-gluon vertices 
exhibit ideal particle-like poles [a = 1 in Eq. 



(2.24)]. Since these poles are an essential element 
of the solution procedure then, in the absence of a 
physically sensible interpretation or explanation 
of them, one could simply reject this result. 
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Alternatively, one can suppose that Eq. (5.4) 
is more robust than the procedure employed to 
motivate it and explore the phenomenological con- 
sequences of the conjecture [34]: 



(5.5) 



where 5 is a dynamically generated mass scale. 
Following this approach it was found that if there 
are no particle-like singularities in the quark-glu- 
on vertex, ry(g,p), then Vs{k^) is unable to con- 
fine quarks [32,46] and b must be fine-tuned to 
very small values: 



& < 6^ ~ A 



QCD : 



(5.6) 



if DCSB is to occur [32,45-48]. It is therefore ap- 
parent that Eq. ( |5.4| ) is phenomenologically dif- 
ficult to maintain. [NB. Achieving DCSB by re- 
quiring ^ is indicative of the dynamical eva- 
sion of Eq. (5.4) since 'Ps{k'^) ^ 1 rapidly for 
small values of b.] 

Nevertheless, the hypothesis has been explor- 
ed in studies [49] of the dressed-gluon 2-point 
function using numerical simulations of lattice- 
QCD. 7^(fc^ — 0) is necessarily finite in simula- 
tions on a finite lattice because of the inherent 
infrared cutoff. Thus one can only truly deter- 
mine V{k^ ~ 0) by considering the behaviour of 
the numerical result in both the countable limit 
of infinitely many lattice sites and the continuum 
limit. 

The form Vs{k^) does not provide as good 
a fit to the lattice data as an alternative form, 
which in the countable limit is 



Vdk^) 



NP + Zk"^ [k^a"^) 



2 yi ' 



(5.7) 



50Aq(~,j-,, where 1/a 



< fc2 < o.6/a2 
2.0 GeV is the inverse lattice spacing, Z « 0.1, 
r] « 0.53, and M « 0.16 GeV. This takes the 
maximum value 

7'L(fc2 = 21A2jcD) = 13.6 (5.8) 

and corresponds to a less extreme alternative to 
Eq. ( |5.4| ), which I shall characterise as 

V{k^ = 0) = , max {Vik^)) < 0(10) . (5.9) 

The feature 7'(fc^ = 0) = is critically depen- 
dent on whether or not M is nonzero. It ap- 
pears to be nonzero in the countable limit but, 



as emphasised in Ref. [49], the behaviour of M 
(and rj) in the continuum limit is unknown. [NB. 
All existing lattice-QCD simulations of the gluon 
propagator; e.g., Ref. [50], yield fitted forms that 



lie in the class specified by Eq. (5.9). A dressed- 
gluon propagator satisfying Eq. (5.4) automati- 



cally satisfies Eq. (p. 9]).] 

The phenomenological implications of Eq. 



(5.7) can be explored using the methods of Ref. 
[32] . A preliminary estimate follows from observ- 
ing that ^'^(fc^) is approximately equivalent to 
Vsik"^) if one identifies 



bL \jMla = 0.57 GeV. 



(5.10) 



Hence one expects that Eq. (5.7) does not gen- 
erate DCSB nor confine quarks. [A value of 6 « 
0.4GeV> be in Vs{k^) provides the best fit to 
the lattice data and this supports the same con- 
clusion.] In order to quantitatively verify this 
conclusion I note that: it is the combination 



g'vik^Mk^ 



(5.11) 



that appears in Eq. (2.1) and is not deter- 
mined in Ref. [49]; and one must extrapolate 
■pL(fc^) outside the fitted domain. Both of these 
requirements are fulfilled if: 1) one assumes that 
a one-loop perturbative analysis is reliable for 



25 Aqqq; and 2) employs 



a'Viik^) := 



,^L(fc2), k^<kl 



gHk'), 



k^ > kl 



(5.12) 



with g^ikF') the one-loop running coupling, re- 
quiring that A;(fc^) := Vi{k'^)/k'^ and its first 
derivative be continuous at k'^-^. This procedure 
yields Ai(fc^) in Fig.|| with 



9r, 



0.65, kt, 



30A|cD- 



(5.13) 



It is now straightforward to solve Eq. (2.1) 
with a variety of Ansdtze for the quark-gluon ver- 
tex that do not exhibit particle-like singularities. 
This was pursued in Ref. [48], using the meth- 
ods described in Sec. ^ and renormalising at the 
momentum cutoff, Auv ~ ICAqcd, for simplic- 
ity, since the p^-evolution of A{p^) and B{p^) 
beyond that point is completely determined by 
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■ lattice gluon model, eq. (9) 
one-loop ACk ) 
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Figure 5: A;(fc2) := Vi[k^)lk^ from Eq. ( jsl^ ). 
Viik^) is Eq. (5.7) in the infrared and extrapolates 



this lattice model outside the domain accessible in 
the simulation [49,50]. (Adapted from Ref. [48].) 

g^{k^). That study employed the bare vertex 
^p.iP,q) 7^; the Ansatz[l4]: 

tT^'^ip, q) zS^(p2, g2) 7, + (p + q), (5.14) 
X [i^7•(p + g)AA(p^g') + As(p^<72)] , 

where 



SA(p^g2) 

Aa(p',9') 
Ab(p2,(z2) 



= [A{p^) + A{q^)]l2, 

= [A{p^) - A{q^)]/[p' - q^] 

^[B{p^)^B{q^)]/[p^ -q'^ 



(5.15) 



and an augmented form [8] 



rCP(p;q):=rBC(p,5) + r6(p,g), (5.16) 



(5.17) 



A{p^) - A{q^) 
2d{p, q) 



with 



d{p, q) 



[p' 



-,212 



+ [M(p2)2+ 1/(92) 



2\2]2 



p2 ^ q2 



(5.18) 

each of which allows the quark DSE to be solved 
in isolation; i.e., without coupling to other DSEs. 



In all cases the lattice result, Eq. (5.12), yields 
M(p2) = in the chiral limit; i.e., no DCSB . 
The absence of DCSB means it is straightfor- 



ward to decide whether ( 5.12 ) generates confine- 
ment. Following Sec. H, confinement is manifest 



Figure 6: Z{p^) obtained as the solution of Eq. ( 2.1) 
using Eq. ( ^.12[ ) with; Eq. (^.14]) - solid line; Eq. 
Aq) - dashed line; and 



7m 



dotted line. 



That Eq. (5.7) does not confine quarks is manifest in 
the result: Z{j)^ = 0) 7^ 0, which is independent of 
the vertex Ansatz. (Adapted from Ref. [48].) 



if, for p2 ~ 0, 



Zip") cc {p^r, a > 1, 



(5.19) 



in which case the dressed-quark propagator does 
not have a Lehmann representation; i.e., violates 
the axiom of reflection positivity. The solution 
obtained with the vertex Ansdtze introduced 
above is depicted in Fig. ^: the behaviour is 
qualitatively equivalent in each case and demon- 



strates explicitly that Eq. (5.7) does not generate 
confinement. 

These failures: the absence of both DCSB 
and confinement, confirm the preliminary hypoth- 
esis based on the correspondence with Vs via an 
effective value of b, Eq. (5.10). [The same con- 



clusion applies to the result in Ref. [50] , which is 
is pointwise smaller in magnitude (< 1/3) than 
Eq. (p??!) on the entire fitted domain.] 



Equation (2.1) was also solved for Z{p'^) us- 



ing 



Piik^) := (1 <j e-'^'/^QCD^ 7?;(fc2) (5 20) 

where ^ is a variable "strength" parameter. In- 
creasing <: moves the peak in Vi^k"^) toward fc^ = 
and increases its height, thereby making it in- 
creasingly like the model of Eq. ( |3.9|) . The form 
of Z{p'^) is qualitatively unchanged and hence 
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there is no signal for the onset of confinement 
until <j > 300. At = 300 the maximum value is 



^/(fc2 = 0.98A|cD) = 210 



(5.21) 



and 



(0.98 A? 



w)/^'(30A^cd) = 16.0, 
of. Vi{2lKlc^)/Vi{iQKlcB) = 1-0 • 



(5.22) 

The model of Eq. (3.£) exhibits a peak at k"^ = 



3.7Aq(;;]3 and the value of the ratio introduced 
above is 44, neglecting only for the purpose of 
this comparison the purely long-range, (5''(fc)-part 
of that interaction. A comparison of 



QCD 



) « 89 



(5.23) 



with the critical coupling of 5^ ~ 11 in Refs. [24] 
shows that such large values of <; ensure DCSB. 

To recapitulate: in this section we have seen 
that in the absence of particle- like singularities in 
the dressed-quark-gluon vertex, extant propos- 
als for the dressed-gluon propagator that satisfy 



Eq. ( 5^ ) neither confine quarks nor break chiral 
symmetry dynamically. This class includes all 
existing estimates of V{k^) via numerical simu- 
lations. 

5.1 Critical Interaction Tension 

Furthermore the calculations described in this 
section reaffirm the long-standing DSE result [24] 
that the existence of DCSB in QCD places a con- 
straint on the minimum phenomenologically ac- 
ceptable value of the "interaction tension:" 

1 r'^pQcu 



/ dfc2fc2e(fc2) (5 24) 



47r ./A 

^^QCD 



The lattice result, Eq. ( 5.12 ), yields 



~0.35GeV2^6A|cD, 



(5.26) 



while from Ref. [32] I estimate the critical value 
for DCSB to be: 



> - 0.5 GeV^ 



(5.27) 



with the exact value depending a little [~ 20%] 
on model details. This is the critical value, which 
only supports an incipient quark condensate; i.e., 
(qq) infinitesimally-greater-than zero. For agree- 
ment with its phenomcnological value one requi- 
res [18]: 



phcn. 



4GeV^ 



70A|cd; 



(5.28) 



i.e., on order-of-magnitude larger. 

One additional observation: following Refs. 
[6,51], an estimate of the gluon condensate can 
be obtained from 



{aGG) := 



47r3 



^pQCD 



dfc2fc4e(fc2)_ (5 29) 



QCD 



The lattice model, Eq. ( ^l2[ ), yields 0.09 GeV; 
the critical value for DCSB is w 0.12 GeV; and 
the model of Ref. [18] yields 0.42 GeV''. Con- 
temporary analyses of QCD Sum Rules employ 
a value [27]: 0.095 ± 0.05 GeV^. However, this 
cannot be directly compared with the calculated 
values reported here because of the ambiguity 
in subtracting the "perturbative contribution" to 



the gluon correlator. Equations (5.24) - (5.2£), 
with a subtraction that is constant below fc^ = 
■^pQCD) employ just one of many possibilities. 
Nevertheless the calculated results indicate the 
relevant scales. 



e(fc^) := g'Aik') ~ g'A{A') ■ (5-25) 6. Gluon DSE 



^pQCD 

The lower cutoff is Aqcd because string-breaking 
becomes effective for lengths > 1 /Aqcd ^ 1 fm, 
while the upper cutoff, ApQco = 2 GeV, marks 
the boundary above which the interaction is cal- 
culable in perturbation theory. [NB. If a mecha- 
nism can be found by which singularities in the 
dressed-quark-gluon vertex become phenomeno- 
logically tenable then this constraint persists with 
the only change being g^A(fc^) G{k'^)/k'^; i.e., 
the effective interaction.! 



Plainly now it is a model-independent result that 
DCSB requires a significant infrared enhancement 
in the kernel of QCD's gap equation. Such an 
enhancement can also yield confinement, by en- 
suring that coloured n-point Schwinger functions 
violate reflection positivity, as discussed in Sec. 
H The question that naturally arises is: "Where 
does this enhancement come from?" Some guid- 
ance may be sought in studies of the DSE sat- 
isfied by the dressed-gluon propagator, which is 
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depicted in Fig. ^. However, as I now describe, 
these studies are inconclusive. 

Early analyses [52] used the ghost-free axial 
gauge: n ■ A°- = > 1^ in which case the 

second equation in Fig. |4| is absent and two inde- 
pendent scalar functions: Fi, F2, are required to 
fully specify the dressed-gluon propagator, cf. the 



covariant gauge expression in Eq. ( 2.17 ), which 
requires only one function. In the absence of in- 
teractions: = -l/P, F2{P) = 0. These 
studies employed an Ansatz for the three-gluon 
vertex that doesn't possess a particle-like singu- 
larity and neglected the coupling to the quark 
DSE. They also assumed F2 = 0, even nonper- 
turbatively, and ignored it in solving the DSE. 
The analysis then yielded 



Flip) ^ ^; 



(6.1) 



i.e., a marked infrared enhancement that can yield 
an area law for the Wilson loop [53] and hence 
confinement, and DCSB as described above, with- 
out fine-tuning. 

This effect is driven by the gluon vacuum po- 
larisation, diagram three in the first line of Fig. ^. 
A similar result was obtained in Ref. [54]. How- 
ever, a possible flaw in these analyses was identi- 
fied in Ref. [55] , which argued from properties of 
the spectral density in ghost-free gauges that F2 
cannot be zero but acts to cancel the enhance- 
ment in Fi. [NB. Retaining F2 in the analysis 
yields a coupled system of equations for the gluon 
propagator that is at least as complicated as that 
obtained in covariant gauges, which perhaps out- 
weighs the apparent benefit of eliminating ghost 
fields in the first place.] 

There have also been analyses of the gluon 
DSE using Landau gauge and those of Refs. 
[56,57] are unanimous in arriving at the covari- 
ant gauge analogue of Eq. (^), again driven by 
the gluon vacuum polarisation diagram. In these 
studies Ansdtze were used for the dressed-three- 
gluon vertex, all of which were free of particle-like 
singularities. However, these studies too have 
weaknesses: based on an anticipated dominance 
of the gluon-vacuum polarisation, truncations 
were implemented so that only the third and fifth 
diagrams on the r.h.s. of the first equation in 
Fig. ^ were retained. In covariant gauges there is 



a priori no reason to neglect the ghost loop con- 
tribution, diagram six, although perturbatively 
its contribution was estimated to be small [57]. 

As described on page the Landau gauge 
studies of Refs. [34] yield a qualitatively differ- 



ent result: Eq. (5.5), but the question of how the 
particle-like singularities in the associated vertex 
Ansdtze can be made consistent with the absence 
of coloured bound states in the strong interaction 
spectrum is currently unanswered. Nonetheless 



proponents of the result in Eq. (5^) claim sup- 
port from studies [58,59] of "complete" gauge 
fixing; i.e., in the outcome of attempts to con- 
struct a Fadde'ev-Popov-like determinant that 
eliminates Gribov copies or ensures that the func- 
tional integration domain for the gauge field is re- 
stricted to a subspace without them. Fixing a so- 
called "minimal Landau gauge," which enforces 
a constraint of integrating only over gauge field 
configurations inside the Gribov horizon; i.e., on 
the simplest domain for which the Faddeev-Popov 
operator is invertible, the dressed-gluon 2-point 
function is shown to vanish at fc^ — 0. However, 
the approach advocated in Refs. [34] makes no 
use of the additional ghost-like fields necessary 
to restrict the integration domain. [NB. Hith- 
erto the quantitative effect of Gribov copies in 
nonperturbative studies remains unknown.] 

Recently the direct approach to solving the 
Landau gauge gluon DSE, pioneered in Refs. 
[56,57], has been revived by two groups: A, Refs. 
[44,60,61]; and B, Refs. [62-65], with the signifi- 
cant new feature that nonperturbative effects in 
the ghost sector are admitted; i.e., a nonpertur- 
bative solution of the DSE for the ghost propa- 
gator is sought in the form 



cab 



fc2 



(6.2) 



[Without interactions, ■cu{k'^) = 1.] 

These studies analyse a truncated gluon- 
ghost DSE system, retaining only the third and 
sixth loop diagrams in the first equation of Fig. ^, 
and also the second equation. Superficially this 
is the same complex of equations as studied in 
Refs. [34]. However, the procedure for solving it 
is different, arguably less systematic but also less 
restrictive. 
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The difference between the groups is that 
A employ Ansdtze for the dressed-ghost-gluon 
and dressed-three-gluon vertices constructed so 
as to satisfy the relevant Slavnov- Taylor identi- 
ties while B simply use the bare vertices. Never- 
theless they agree in the conclusion that in this 
truncation the infrared behaviour of the gluon 
DSE's solution is determined by the ghost loop 
alone: it overwhelms the gluon vacuum polarisa- 
tion contribution. 

That is emphasised in Ref. [64], which elim- 
inates every loop diagram in truncating the first 
equation of Fig. ^ except the ghost loop and still 
recovers the behaviour of Ref. [65]. That be- 
haviour is 



1 



(6.3) 



for k'^ < ^^QD, with 0.8 < k < 1. Exact eval- 
uation of the angular integrals that arise when 
solving the integral equations gives the integer- 
valued upper bound, k = 1 [63]. This corre- 
sponds to a dressed-gluon 2-point function that 
vanishes at fc^ = 0, although the suppression is 
very sudden with the propagator not peaking un- 
til fc^ w Aqqq, where 



{d{e)ie)\ 



fe2=A 



QCD 



lOO/A^ 



QCD : 



(6.4) 



i.e., it is very much enhanced over the free prop- 
agator. [See, e.g., Ref. [44], Fig. 12.] k = 1 
also yields a dressed-ghost propagator that ex- 
hibits a dipole enhancement analogous to that of 



Eq. ( |6T| ). 

A [renormalisation group invariant] running 
strong-coupling consistent with the truncations 
that yield these solutions is: 



(6.5) 



and its value at fc^ = is fixed by the numerical 
solutions: 



(6.6) 





A 


B 


a{e = 0) 


9.5 


~ 4 or 12 



[NB. Group A approximates the angular inte- 
grals and uses vertex Ansdtze. Group B uses 
bare vertices and in Ref. [62] approximates the 
angular integrals to obtain q;(0) ~ 12, while in 



Ref. [63] the integrals are evaluated exactly, which 
yields a(0) = |7r « 4.2.] 

The qualitative feature common to both 
groups is that the Grassmannian ghost loops act 
to suppress the dressed-gluon propagator in the 
infrared. That may also be said of Refs. [34]. 
[Indications that the quark loop, diagram seven 
in Fig. 0, acts to oppose an enhancement of the 



type in Eq. (6.1) may here, with hindsight, be 
viewed as suggestive.] 

One aspect of ghost fields is that they enter 
because of gauge fixing via the Fadde'ev-Popov 
determinant. Hence, while none of the groups 
introduce the additional Fadde'ev-Popov contri- 
butions advocated in Refs. [58,59], they neverthe- 
less do admit ghost contributions, and in their so- 
lution the number of ghost fields does not have a 
qualitative impact. Reference [59] also obtains a 
dressed-propagator for the Fadde'ev-Popov fields 
with a fc^ = dipole singularity. It contributes 
to the action via the term employed to restrict 
the gauge field integration domain, in which ca- 
pacity the dipole singularity can plausibly drive 
an area law for Wilson loops. 

Schwinger functions arc the primary object 
of study in numerical simulations of lattice-QCD 
and Refs. [49,50] report contemporary estimates 
of the lattice Landau gauge dressed-gluon 2-point 
function. As we saw in Sec. they are consis- 
tent with a finite although not necessarily van- 
ishing value of (i(fc^ = 0). However, simulations 
of the dressed-ghost 2-point function find no evi- 
dence of a dipole singularity, with the ghost prop- 
agator behaving as if w{k'^) = 1 in the smallest 
momentum bins [67]. [NB. Since the quantita- 
tive results from groups A and B differ and ex- 
hibit marked sensitivity to details of the numer- 
ical analysis, any agreement between the DSE 
results for Tu{k'^) or d{k'^) and the lattice data on 
some subdomain can be regarded as fortuitous.] 

The behaviour in Eqs. (^^) also entails the 
presence of particle-like singularities in extant 
Ansdtze for the dressed-ghost-gluon, dressed- 
three-gluon and dressed-quark-gluon vertices that 
are consistent with the relevant Slavnov- Taylor 
identities, [k = 1 corresponds to an ideal simple 
pole singularity.] Hence while this behaviour may 
be consistent with the confinement of elementary 
excitations, as currently elucidated it also pre- 
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diets the existenee of eoloured bound states in 
the strong interaction spectrum [see page |l^ . 

Furthermore, while it does yield a running 
strong-coupling with a(fc^ =0) ^ 1) that, as 
we saw in Sec. ||, makes DCSB dependent on 
fine tuning. To make this plain, the maximum 
value of 12 in Eq. (U) yields cr^ w 0.7 GeV^, 
which is just above but falls far short of the 
value required to produce the physical value of 
the quark condensate. In fact, the quark conden- 
sate is only ~ 5% of the phenomenologically re- 



quired value in Eq. ( 3.3C ) [66] . Notwithstanding 
these remarks, the studies of Refs. [44,60,61] and 
subsequently Refs. [62-65] are laudable. They 
have focused attention on a previously unsus- 
pected qualitative sensitivity to truncations in 
the gauge sector. 

To recapitulate. It is clear from Sec. || that 
DCSB requires the effective interaction in the 
quark DSE to be strongly enhanced at ~ 
^QCD- [Remember too that in Sec. ||, Figs. 
we saw that modern lattice simulations confirm 
the pattern of behaviour exhibited by quark DSE 
solutions obtained with such an enhanced inter- 
action.] Studies of QCD's gauge sector indicate 
that gluon-gluon and/or giuon-ghost dynamics 
can generate such an enhancement. However, 
the qualitative nature of the mechanism and its 
strength remains unclear: is it the gluon vac- 
uum polarisation or that of the ghost that is the 
driving force? It is a contemporary challenge to 
explore and understand this. 



7. Bethe-Salpeter Equation 

Hitherto I have focussed on the not-directly-ob- 
servable elementary excitations in QCD. They 
make themselves manifest in hadrons and their 
properties. In quantum field theory, two and 
three-body bound states are, respectively, des- 
cribed by the Poincare-covariant Bethe-Salpeter 
and Fadde'ev equations. Solving an equation 
of this type yields the bound states mass and 
also an amplitude that describes the bound state 
constraints on the constituents' momenta. This 
amplitude is a valuable intuitive guide and, in 
cases where a simple quark model analogue of 
the bound state exists, the amplitude incorpo- 



rates and extends the information present in that 
analogue's quantum mechanical wave function. 

In quantum field theory, as in classical me- 
chanics, the interacting two-body problem is 
much simpler than that of three such bodies, and 
I illustrate the bound state application of DSEs 
by considering mesons. The renormalised homo- 
geneous BSE for a bound state of a dressed-quark 
and dressed-antiquark with total momentum P is 



[ri/(fc;P)]* 



[x{q-,P)]rsKZ{q,k;P)., 



(7.1) 

XH{q; P) ■■= Siq+) THiq; P) S{q+) , (7.2) 

with: TH{k;P) the Bethe-Salpeter amplitude, 
where H specifies the flavour structure of the 
meson; S{p) := dia,g[Su{p), Sd{p), Ss{p)]; q+ = 
q + rjpP, q^ ~ q ~ (1 — •qp)P\ and r, . . . ,u rep- 
resent colour-, Dirac- and flavour-matrix indices. 
[rjp € [0, 1] is the momentum partitioning pa- 
rameter. It appears in Poincare covariant treat- 
ments because, in general, the definition of the 
relative momentum is arbitrary. Physical observ- 
ables, such as the mass, must be independent of 
rjp but that is only possible if the Bethe-Salpeter 
amplitude depends on it. 77P = 1/2 for charge- 
conjugation eigenstates.] 

The general form of the Bethe-Salpeter am- 
plitude for pseudoscalar mesons is 



VH{k;P) = 



(7.3) 



iEH{k-P)+^-PFH{k;P) 



-f-kk-P Gnik; P) + a^,, k^P, Hnik; P) 



where for bound states of constituents with equal 
current-quark masses, the scalar functions E, F, 
G and H are even under k ■ P —k ■ P. [NB. 
Since the homogeneous BSE is an eigenvalue 
problem, EH{k;P) = Enik^ , k ■ P\P^), etc.; i.e., 
P^ is not a variable, instead it labels the solu- 
tion.] 



In Eq. (7.1), Kl^{q, k; P) is the renormalised, 
fully-amputated quark-antiquark scattering ker- 



nel, which also appears implicitly in Eq. (2.1) 
because it is the kernel of the inhomogeneous 
DSE satisfied by T^{q;p). Kl%q,k;P) is a 4- 
point Schwinger function, obtained as the sum 
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of a countable infinity of skeleton diagrams. It 
is two-particle-irrcducible, with respect to the 
quark- antiquark pair of lines and does not con- 
tain quark-antiquark to single gauge-boson anni- 
hilation diagrams, such as would describe the lep- 
tonic decay of a pseudoscalar meson. (A connec- 
tion between the fully-amputated quark-antiquark 
scattering amplitude: M = K -\- K{SS)K -I- . . . , 
and the Wilson loop is discussed in Ref. [68].) 

The complexity of Kl^{q, k;P) is one reason 
why quantitative studies of the quark DSE cur- 
rently employ Ansdtze for D^i,[k) and r^(fc,p). 
However, as illustrated by Ref. [25], this com- 
plexity does not prevent one from analysing as- 
pects of QCD in a model independent manner 
and proving general results that provide useful 
constraints on model studies. [NB. References 
[25,18,69-71] provide a pedagogical guide to the 
rigorous and practical application of BSEs to the 
light-quark sector of QCD, and should be par- 
ticularly helpful to practitioners for whom their 
efficacy has not hitherto been apparent.] 

K has a skeleton expansion in terms of the 
elementary, dressed-particle Schwinger functions; 
e.g., the dressed-quark and -gluon propagators. 
The first two orders in one systematic expan- 
sion [3] are depicted in Figs. 0, ||. This partic- 
ular expansion, in concert with its analogue for 
the kernel in the quark DSE, provides a means 
of constructing a kernel that, order- by-order in 
the number of vertices, ensures the preservation 
of vector and axial-vector Ward-Takahashi iden- 
tities. This is particularly important in QCD 
where the Goldstone boson nature of the pion 
must be understood as a consequence of its inter- 
nal structure, and without fine-tuning; i.e., the 
masslessness of the pion in the chiral limit can- 
not arise as the result of carefully varying param- 
eters in a putative potential. The Goldstone bo- 
son character of the pion is easily understood via 
the interrelation between the pseudoscalar BSE 
and quark DSE [3,4,25]. 

Analysing the flavour-nonsinglet inhomoge- 
neous axial- vector BSE and its pseudoscalar ana- 
logue one obtains, via the axial-vector Ward-Ta- 
kahashi identity, a mass formula, exact in QCD, 
for pseudoscalar mesons [4,25,72]: 




Figure 7: Leading [Ladder] contribution to the sys- 
tematic expansion of the quark-antiquark scattering 
kernel introduced in Ref. [3]. In this expansion, the 
propagators are dressed but the vertices are bare. 
(Adapted from Ref. [3].) 



= trF[MjT^, (T^)'}], where is a ma- 
trix that describes the flavour content of the me- 
son; e.g., = ^{Xp+iXp), (•)* denotes matrix 
transpose, and 



fHPfi = Zi tr 



T 



' H 



(7.5) 



= -iV2Zitr {T") 75 XH (fc; P) 



(7.6) 



Importantly, the expression in Eq. ( [7.4[ ) is valid 
for all values of the current-quark mass. [NB. In 
this section the normalisation is such that /^r = 
131 MeV.] 



Equation (7.5) gives the formula for the re- 
sidue of the meson's pole in the axial- vector ver- 
tex, which completely determines the strong in- 
teraction contribution to its leptonic decay. The 
factor of Z2 on the r.h.s. is just that necessary 
to ensure that fn is independent of the renor- 
malisation point, regularisation mass-scale and 
gauge parameter; i.e., to ensure that fn is truly 
a physical observable. Its intuitive character is 
also plain: it is the gauge-invariant projection of 
the meson's Bethe-Salpeter wave function at the 
origin. 



(7.4) 



Equation (7^) is the expression for the resi- 
due of the meson's pole in the pseudoscalar ver- 
tex. In this case again the factor Z4 on the r.h.s. 
depends on the gauge parameter, the regulari- 
sation mass-scale and the renormalisation point. 
This dependence is exactly that required to en- 
sure that: 1) rn is finite in the limit A — s- 00; 2) 
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Figure 8: Next-to-leading order contribution in the truncation scheme of Ref. [3]. (Adapted from Ref. [3].) 



ru is gauge-parameter independent; and 3) the 
renormalisation point dependence of th is just 
such as to ensure that the r.h.s. of Eq. ( |7.4D is 
renormaUsation point independent. 

These formulae exemplify the manner in 
which gauge invariant results are obtained from 
the gauge covariant Schwinger functions the DSEs 
yield. 

It is straightforward to show that in the chi- 
ral limit, defined in Eq. (2.15), Eq. (^) yields 



1.0 



fO 



(7.7) 



where (qq)^ is the vacuum quark condensate, Eq. 
( 3.26| ), and /" is the pseudoscalar meson decay 
constant in the chiral limit. Hence, as a corollary 



of Eq. (7.4) one recovers [25] what is commonly 
called the "Gell-Mann-Oakes-Renner" relation: 



f^mjj^im^ (qq)° 



(7.8) 



Another particularly important result is that 
the axial- vector Ward-Takahashi identity also con- 
strains the chiral limit behaviour of the scalar 



functions in Eq. (7.3). This is a pointwise mani- 
festation of Goldstone's theorem with [25]: 

^/Oi?4fc;0) = i?(fc2), (7.9) 

i^H(fc;0) + y2/X(fc;0) = A(fc2), (7.10) 

Gfl(fc;0) + V2/°G,(fe;0) = 2A'(fc2),(7.11) 

i?fl(A:;0)+V2/0F^(fc;0) = 0. (7.12) 

In these identities (•)^ denotes solution functions 
obtained in the chiral limit and [-)r are functions 
appearing in the pole-free part of the inhomoge- 
neous axial- vector vertex. 



The formula in Eq. (7.£ ) is a quark-level Gold- 
berger-Treiman relation, exact in QCD, while the 
next two identities indicate that pseudoscalar me- 
sons necessarily have a pseudovector component. 



0.6 



0.4 



0.2 



0.0 



0.0 



0.2 



0.4 
q' (GeV") 



0.6 



Figure 9: Small-g^ pion form factor calculated in 
Ref. [73] and compared with data from Refs. [75] 
(circles) and [76] (crosses). The solid line is the re- 
sult obtained including the pseudovector components 
while the dashed line was obtained without them 
[74]. Clearly they contribute little on this domain. 
(Adapted from Ref. [73].) 



This is crucial because it is these amplitudes that 
dominate the mesons' electromagnetic form fac- 
tors at large meson energy and yield [73] 



q^Fniq^) — constant at large— 



(7.13) 



up to logarithmic corrections, making a connec- 
tion with the result of perturbative QCD analy- 
ses. If the pseudovector amplitudes are ignored 
one finds instead [74] that (g^)^-Ff/(g^) = con- 
stant at large-q^. These results are illustrated in 
Figs. ||, pO|. They are, of course, contingent 



upon Eq. ( 3.25| ), which describes the correct ul- 
traviolet behaviour of the mass function in QCD. 
Any model that employs/predicts a mass func- 
tion that falls faster than a single power of 
will be in conflict with perturbative analyses of 
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6.0 



Figure 10: Calculated [73] form factor compared 
with the largest (f' data available: diamonds - Ref. 
[75]; and circles - Ref. [77]. The solid lines are the 
results obtained when the pseudovector components 
of the pion Bethe-Salpeter amplitude are included 
(two limiting models were employed in the calcula- 
tion), the dashed-line when they are neglected [74]. 
On this domain the difference is plain, with only the 



solid lines exhibiting Eq. (7.13). (Adapted from Ref. 
[73].) 



Figure 11: Mass-dependence of the meson leptonic 
decay constant. Experimental values of /tt.k, filled 
circles; lattice estimates [79] of fo = 0.20 ±0.03 GeV 
and /s = 0.17 ± 0.035 GeV, open circles; decay con- 
stants calculated in Ref. [22] with an estimated 10% 
theoretical error. The dashed line is an eye-guiding 
fit to the experimental and lattice estimates, which 
exhibits the large-rrii/ limit of Eq. (7.17). The dot- 
ted line is the large-m_r/ limit of this fit. (Adapted 
from Ref. [22].) 



the form factors. [NB. No bound state amplitude 
can fall slower than l/p^] 

7.1 Heavy-quark Limit 

After this focus on the light-quark sector I now 



return to Eqs. ( [7.4[ )-(7.6) and note that they also 
have two important corollaries valid in the limit 
of large current-quark mass. To elucidate them 
one rewrites the total momentum 



P^L = mHVfj^ = [Mq + Eh) 



(7.14) 



where: mu is the heavy-hadron's mass; Mq is a 
constituent-heavy-quark mass, Mq m Mq « toq 
for heavy-quarks, as illustrated in Fig. Q and Eh 
is the "binding energy." 

Now the dressed propagator for the heavy- 
quark in the heavy-meson can be written [30] 



SQ{k + P) 




2 k-v- Eh 



where k is the momentum of the lighter con- 
stituent, and the canonically normalised Bethe- 
Salpeter amplitude can be expressed as 



TH{k;P) 



/7^T<r{k-P) 



(7.16) 



where T'^°°{k]P) is pointwise-finite in the limit 
niH — > oo. It is obvious that in the calcula- 
tion of observables the meson's Bethe-Salpeter 
amplitude will limit the range oi\k\ so that Eq. 



(7.15) will be a good approximation if-and-only- 
if both the momentum space width of the am- 
plitude and the binding energy are significantly 
less-than Mq. 

Using Eqs. ( [7T5| ) and ( |7l6| ) in Eq. ( fT^Sj ) 
yields 

fH = ^. (7.17) 



where c^^j is a calculable and finite constant that 
depends only on r^°°(fc; P) and the dressed-light- 



quark propagator. In Eq. (7.17) the DSEs repro- 
duce a well-known general consequence of heavy- 
quark symmetry. Since /tt < fx is an experimen- 
tal fact neither this formula nor the limits that 
lead to it are valid for light mesons. The actual 
situation, as determined in Ref. [22], is depicted 
in Fig. |l^. Clearly, finite current-quark mass cor- 
rections are significant for the c-quark. In fact, 
such corrections can be < 30% in 6 c transi- 
tions and as much as a factor of two in c — s- d 
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Figure 12: Mass-dependence of meson masses. Ex- 
perimental value of the masses and estimates of niq, 
are taken from Ref. [17]. The straight line, drawn 
through the D- and B-meson masses, depicts the lin- 
ear trajectory predicted in Eq. (7.19). (Adapted from 
Ref. [78].) 



transitions [22]. 

Applying the same analysis to Eq. ( |7.6| ) gives 

TH = c'h v^'^h , (7-18) 

where again cj^ is a calculable and finite con- 
stant. Using this and Eq. ( [7.17 ) in Eq. (7.4) one 
finds [22,78] that in the heavy-quark limit 



UlH 



(7.19) 



Thus the quadratic trajectory of Eq. (7.8) evolves 
into a linear trajectory when the current-quark 
mass becomes large; i.e., the mass of a heavy- 
meson rises linearly with the mass of its heaviest 



constituent. This is illustrated in Fig. 12, which 
indicates that a linear trajectory may describe all 
but the TT-meson. That was anticipated in Ref. 
[4] and later clarified in an explicit calculation 
[21] where the onset of the linear trajectory was 
observed at 2~3-times the s-quark current-mass. 

Explicit calculations, of course, require that 
the form of the scattering kernel be specified. 
The dressed-ladder form 

Kl^{q,k;P)= (7.20) 



in concert with Eq. ( |3.9| ) in the quark DSE, en- 
sures that the vector and axial-vector Ward-Ta- 



kahashi identities are satisfied. Preserving these 
identities is important because it guarantees cur- 
rent conservation in electromagnetic processes 
and an explicit realisation of the consequences 
of DCSB as elucidated herein. The choice in Eq. 



( |7.20| ) is that yielded by the truncation scheme 
of Ref. [3] and clearly maintains the ultraviolet 



behaviour of QCD, Eq. ( |3.5D . 

With an explicit form for the kernel, Eq. 



(7.1) is an eigenvalue problem and it has solu- 
tions only for particular, separated values of P^. 
The eigenvector associated with each eigenvalue: 
Tnik; P), the Bethe-Salpeter amphtude, is a one- 
particle-irreducible, fully-amputated quark- 
meson vertex. In the flavour non-singlet pseu- 
doscalar channels the solutions having the lowest 
eigenvalues correspond to the tt- and A'-mesons, 
while in the vector channels they correspond to 
the LU-, p- and 0-mesons. 

To calculate meson masses, one first solves 
Eqs. ( 3.14| ) and ( 3.15| ) for the renormalised dres- 
sed-quark propagator. This numerical solution 
for S{p) is then used in the BSE [obtained using 



Eqs. ([73D, (Q, ( [7.20D ], which for pseudoscalar 
mesons is a coupled set of four homogeneous equa- 
tions, one set for each meson. Solving the equa- 
tions is a challenging numerical exercise, requir- 
ing careful attention to detail, and two comple- 
mentary methods were both used in Refs. [4,18]. 
While the numerical methods were identical, the 
authors of Ref. [18] used a simplified version of 
the effective interaction: 



-47r ■ 



7m TT 



(7.21) 



1 



P/A^ 



QCD 



and varied the single parameter D along with the 
current-quark masses: rhu — rhd and rhs, in or- 
der to reproduce the observed values of m.^ , rrifc 
and /tt listed on page ||. The many other calcu- 
lated results in Refs. [18,69,70] are predictions in 
the sense that they are unconstrained. A partic- 
ular feature of these studies is that the Poincare 
covariant four-dimensional BSE is solved directly 
eschewing the commonly used artefice of a three- 
dimensional reduction, which introduces spuri- 
ous effects when imposing compatibility with 
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Goldstone's theorem and can also lead to the mis- 
interpretation of a model's parameters [80]. 

8. Epilogue 

I have focused on some robust qualitative as- 
pects of Dyson-Schwinger equation [DSE] stud- 
ies. A compelling complement is Ref. [71], which 
describes in detail the direct application of the 



single-parameter DSE- model of Eq. (7.21) to the 
spectrum and dynamical properties of light me- 
sons. Additional important aspects are described 
in Ref. [81]: in particular the use of vector meson 
polarisation observables to probe the long-range 
part of the quark-quark interaction, which high- 
lights the intuitive character of Bethe-Salpeter 
amplitudes; and also in Refs. [73,74,82] where 
the manner in which the DSEs provide for the 
intrinsic preservation of current-algebra's anoma- 
lies is elucidated and exemplified. Light-baryon 
properties too have been studied, with Ref. [83] 
demonstrating that a quark-diquark simplifica- 
tion [84] of the covariant baryon Fadde'ev equa- 
tion is a good foundation for spectroscopy, and 
Refs. [85] proving it efficacious for a wide range 
of scattering observables. 

In recent years there has also been a renewal 
of interest in DSE applications to QCD at non- 
zero chemical potential and temperature, which 
herein I have not mentioned at all. In closing 
therefore I note that, as will have become clear 
from Sees. |3| and ^, confinement and dynamical 
chiral symmetry breaking are simultaneously ac- 
cessible using the DSEs. This was exploited in 
Ref. [86] wherein a well-constrained DSE-model 
of QCD was used to explore the deconfinement 
and chiral symmetry restoration transitions at 
T 0: in this case the transitions are coincident 
and second order, which is easily understood via 
the causal relation between these two phenom- 
ena and the behaviour of the dressed-quark mass 
function. 

More recently the nonzero temperature and 
density applications have expanded to include 
hadron properties, as illustrated in Refs. [87], 
and have become more refined. A large part of 
Ref. [1] is a chronicle of this. Reference [1] also 
highlights a contemporary challenge in this area: 
the development of a DSE-based quantum trans- 



port theory. That quest is driven by a desire to 
incorporate the treatment of out-of-equilibrium 
aspects of a relativistic heavy ion collision into 
a framework that reliably describes the chemi- 
cal and thermal equilibrium features of hot and 
dense strongly-interacting matter. 
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A. Euclidean Metric Conventions 

For 4- vectors a, b: 



(A.l) 



so that a spacelike vector, Q^, has > 0. The 
Dirac matrices are Hermitian and defined by the 
algebra 

{7^,74 = 2(5^,; (A.2) 



and 



75 := - 71727374 



(A.3) 



so that 

tr [757p7i^7p7o-] = -^^Sfj^i^po- , £1234 1 • (A. 4) 

The Dirac-like representation of these matrices 
is: 



-if 
if 



where the 2x2 Pauli matrices are: 



1 

) - 

1 



-1 



(A.6) 
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